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It is known that when the steady-state of a one-dimensional multi-species system, which evolves 
via a random-sequential updating mechanism, is written in terms of a linear combination of Bernoulli 
shock measures with random walk dynamics, it can be equivalently expressed as a matrix-product 
state. In this case the quadratic algebra of the system always has a two-dimensional matrix represen- 
tation. Our investigations show that this equivalence exists at least for the systems with determin- 
istic sublattice-parallel update. In this paper we consider the totally asymmetric simple exclusion 
process on a finite lattice with open boundaries and sublattice-parallel update as an example. 
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I. INTRODUCTION 

Recently, the steady-state properties of exclusion pro- 
cesses which belong to the class of driven-diffusive sys- 
tems have been under much investigation because of 
their unique physical characteristics such as shocks and 
non-equihbrium phase transitions [H-0. The matrix- 
product approach has turned out to be one of the most 
powerful techniques in determining the steady-states of 
such stochastic non-equilibrium systems which have been 
used to model biological transport and traffic flow. The 
matrix-product approach has been interpreted from dif- 
ferent point of views (a recent review of this approach 
can be found in Q). According to this approach the 
non-equilibrium steady-state weight of any given configu- 
ration of a one-dimensional stochastic system can be con- 
sidered as a matrix element of product of non-commuting 
operators, one for each lattice site, chosen according to 
the state of the lattice site. In order to calculate these 
weights one needs to know a set of algebraic relations 
between these operators. Whether these operators have 
matrix representations is a challenging issue. 

We have recently investigated the relation between the 
dimensionality of the matrix representation of the algebra 
of a given one-dimensional driven-diffusive system with 
nearest-neighbor interactions and the possibility that the 
steady-state of the system in question can be written in 
terms of a linear superposition of product shock mea- 
sures. In this context a shock is defined as a sharp dis- 
continuity in the density profile of particles on the lat- 
tice. In [8| we have shown that for the one-dimensional 
driven-diffusive systems defined on a discrete lattice of 
finite size in which a single product shock measure has 
a simple random walk dynamics under the time evolu- 
tion generated by the stochastic Hamiltonian of the sys- 
tem, the steady-state of the system can easily be writ- 
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ten as a linear combination of these single product shock 
measures. In most cases it is necessary that some con- 
straints on the microscopic reaction rates of the system 
are satisfied. Surprisingly we have seen that at the same 
time the steady-state of the system can be written as a 
matrix-product state with two-dimensional matrix rep- 
resentation for the algebra of the system, provided that 
the same constraints (if they exist) on the microscopic 
reaction rates are satisfied. The matrix representation in 
this case obeys an special structure which contains very 
detailed information about the hopping rates of the shock 
front, as well as the density of particles on the left and 
the right hand sides of the shock front. In present work 
we aim to investigate the same issue; however, this time 
for the systems in a different updating scheme i.e. the 
sublattice-parallel dynamics. 

One of the simplest, yet interesting, driven-diffusive 
model which has been studied widely during recent years 
is the Asymmetric Simple Exclusion Process (ASEP). In 
this model the classical particles enter the system from 
the left boundary of a discrete lattice, diffuse in its bulk 
and leave it from the right boundary with certain rates. 
The derivation of the matrix-product representation from 
the algebraic Bethe ansatz for this model has been stud- 
ied in 0. For this model the equivalence between the 
partition function of the system with random sequential 
dynamics and the partition function of a two-dimensional 
lattice path model of one-transit walks or Dyck paths has 
been studied in 6] . Under the parallel dynamics the par- 
tition function of the ASEP can be expressed as one of 
several equivalent two-dimensional lattice path models 
involving weighted Dyck paths Q. The dynamics of a 
single shock front in the ASEP with a discrete time up- 
dating scheme defined on an infinite lattice has already 
been studied in Q. 

In this paper we will answer the question whether 
the existence of a two-dimensional representation for the 
quadratic algebra of a driven-diffusive system with a dis- 
crete time updating (more specifically sublattice-parallel 
updating) and nearest-neighbor interactions is related 
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to the fact that the steady-state of the system can be 
constructed in terms of a Hnear superposition of prod- 
uct shock measures with simple random walk dynam- 
ics. We will consider the Totally Asymmetric Simple 
Exclusion Process known as TASEP with open bound- 
aries as a simple example. In this model the particles 
only hop towards the right boundary after being injected 
into the lattice from the left boundary. As we will see, 
quite similar to the case of random sequential updat- 
ing scheme studied in it seems that whenever the 
quadratic algebra of the system (since we only consider 
the systems with nearest-neighbors interactions) has a 
two-dimensional matrix representation with a specific 
structure (which will be discussed later) then we can con- 
clude that the steady-state of the system is made up of 
a linear combination of product shock measures with a 
shock front which has simple random walk dynamics and 
vice versa. One of the differences here is that in the 
case of sublattice updating scheme one should define two 
different shocks which behave differently at even or odd 
lattice sites. This has already been observed in [9] for an 
infinite system. 

In the following sections we will first present the math- 
ematical tools and definitions. Then we will study the 
time evolution of two product shock measures defined at 
even and odd sites under the sublattice-parallel update. 
Then we will construct the steady-state of the system in 
terms of a linear combination of these shocks. We will 
bring the quadratic algebra of the system and its two- 
dimensional representation in terms of the shock char- 
acteristics. The conclusion will be presented in the last 
section. 



II. MATHEMATICAL PRELIMINARIES 

Let us start with the definitions and mathematical pre- 
liminaries. Consider a general two-state driven-diffusive 



As can be seen, the particles in the bulk of the lattice 
move only to the right deterministically while obeying the 
exclusion principle. The particles can enter (leave) the 
lattice only from the left (right) boundary with the prob- 
ability a (/3). The dynamics and also the steady-state of 
the TASEP with open boundaries has been proposed and 
studied in 10]. Later its steady-state was studied using 
a matrix formalism in pj| . In the thermodynamic limit 
i.e. L >> 1 one finds that the system has two different 



system with nearest-neighbor interactions and sublattice- 
parallel dynamics in which classical particles move on a 
one-dimensional lattice of length 1L with open bound- 
aries. The bulk dynamics consists of two half time steps. 
In the first half time step even lattice sites and also the 
first and the last lattice sites are updated. From the first 
and the last lattice sites the particles can be injected or 
extracted with certain probabilities. In the second half 
time step only the odd lattice sites are updated. The 
corresponding transfer matrix T consists of two factors 
T = TaTi defined as [H]: 

where T, C and TZ are the matrices for bulk interactions, 
particle input and particle output respectively. The time 
evolution of the probability distribution is governed by 
the following equation: 



T\P[t)) = \P{t + l)). ^ (1) 
In long-time limit the system approaches its steady-state 
and its non-equilibrium probability distribution satisfies 
the following equation: 



T\P*) = \P*). (2) 



As a simple example consider the TASEP which in 
an appropriate basis the above transfer matrix can be 
written as pT| : 



(3) 
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phases: a high-density phase for a > (3 and a low-density 
phase for a < p. A first-order phase transition also takes 
place at the transition line a = f3. 



III. TEMPORAL EVOLUTION OF SHOCKS 

In what follows we study the time evolution of two 
product shock measures using the time evolution equa- 
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tion ([T]). We consider a discrete lattice of length 2L 2k (k = 1, ■ ■ ■ , L) as |/i2fc) and at odd sites 2fc + 1 
and introduce two product shock measures at even sites (/c 0, • • • , L) as |/i2fe+i) according to: 
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We should explain a couple of points here. Firstly, one 
should note that the shock front for |p2fc) hes between the 
lattice sites 2fc— 1 and 2k while the shock front for \fj,2k+i) 
lies between the lattice sites 2k and 2A: + 1. Secondly, 
the shock \h2L+i) indicates aflat distribution of particles 
with densities p"'^'^ and pl'"^" at odd and even lattice 
sites respectively. In this case the shock front can be 



considered to be between the lattice sites 2L and 2L + 1 
which the later is considered as an auxiliary site. Let 
us consider |p2fc) and \p2k+i) as two initial probability 
distributions and investigate their time evolutions using 
([T]) . Suppose that under some constraints on the reaction 
rates the dynamics of these two shocks are given by: 



T\fi2k) = Si\p2k-i) + Sr\p2k+i} + Ss\f^2k) for l<k<L 

T\^J,2k+l) = Slds\p2k) + SrSs\fi2k+2) + <5^|M2/c+3> + Sf\p2k-l) + (Ss + 2SlSr)\^J,2k+l) for 1 < k < L ~ 1 

(6) 

r|pi) = SrSs\fi2) + Sf\p3) + (1 - 5rds - <52)|^i> 

r|p2L + l) SiSs\fi2L) + Sf\p2L-l) + (1 - SlSs - Sf)\p2L + l) 

I 



in which 6s = 1 — St — 6r. 

The first two equations in ^ have already been ob- 
tained in for the ASEP on an infinite lattice; however, 
since the definition of the transfer matrix in this paper 
is based on [ll| , these two time evolution equations have 
exchanged places as compared to the corresponding equa- 
tions (44) and (48) in [9]. For the TASEP with open 
boundaries these equations of motion are valid, provided 
that we have: 

pfd = , pf'^ = 1-/3 

Sr = P , 5i ^ a. 

Note that ([6|) gives a closed set of time evolution equa- 
tions for |/ifc)'s in which A; = 1, • • • , 2i 4- 1. It is also 
interesting to note that in contrast to the ASEP, here 
there is no constraint on the microscopic reaction rates 
(the boundary rates a and (3). Although the particles 
move deterministically towards the right boundary, the 
shock fronts hop both to the left and to the right. This 
is a direct result of the updating scheme. 



IV. STEADY-STATE 

The simplicity of the time evolution equations ([6]) al- 
lows us to construct the steady-state of the system \P*). 
As can be seen, they are similar to the time evolution 
equations for a simple random walker moving on a finite 
lattice with reflecting boundaries; however, the random 
walker behaves differently when it lies at an even or an 
odd lattice site. By considering a linear superposition of 
the shocks as: 

-, 2L+1 

\n = ^ ^>'\t'^) (8) 

k=l 

and requiring that ([2|) should be satisfled one flnds: 

C2fe = ^.(^)''-'forfc = l,--- ,L (9) 
C2fe+i = (^)''forfc = 0,--- ,L. (10) 



4 



The normalization factor or the partition function of the 
system Z can be easily calculated as: 



2lk=i Cfc 



(11) 



Since the steady-state of the system is unique, if one 
calculates the steady-state probability distribution of the 
system using the matrix-product approach, one should 
find the same distribution as (El). 



V. MATRIX PRODUCT APPROACH 

Let us now investigate the steady-state probability dis- 
tribution function of our general two-state model with 
sublattice-parallel dynamics and nearest-neighbor inter- 
actions using the matrix-product approach. According to 
this approach (and in this particular updating scheme) 
the steady-state probability distribution fmiction of the 
system can be written as [11]: 



\V)) (12) 



in which the operators E and D (E and D ) stand for 
the presence of a hole and a particle at odd (even) sites 



respectively. The normalization factor Z is usually called 
the partition function and can easily be written in a grand 
canonical ensemble as: 



i^{W\{E + b)^{E + D)^\V)). 



(13) 



The operators (E, D) and {E, D) besides the vectors \ V)) 
and {{W\ are acting in an auxiliary space. According 
to the standard matrix-product approach by requiring 
that is satisfied one finds that the above mentioned 
operators and vectors should satisfy a quadratic algebra 
given by [lH : 
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Surprisingly, one can see that the following two- 
dimensional matrix representation which can be written 
in terms of the shock hopping rates and the densities of 
the Bernoulli measures at the left and the right hand 
sides of the shock position can generate exactly the same 
probability distribution ([8|: 
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{{W\ = {W^,W2) , \V)) = 
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provided that we have: 



ViW2d = 



-pr°") 



VlW2d : 



(16) 



These relations are nothing but two constraints on the 
parameters vi, V2t wi, 'W2, d and d; therefore, only four of 



these parameters are free. Note that the densities in the 
shock measures and also the shock front hopping rates in 
([6|) should be fixed by the boundaries and the microscopic 
reaction rates; therefore, are not free parameters. 

Let us go back to our simple example. It is shown in 
[m that the TASEP has a quadratic algebra which can 
be written as: 



J 



[E, E] = [D,D]=Q , Ed [E, D], DE = 0, {{W\E{1 ~ a) = {{W\E 
{{W\{aE + D) = {{W\D , (1 - (3)D\V)) = D\V)) , {E + I3D)\V)) = E\V)). 



(17) 
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In the same reference it has been shown that pT)) has 
a two-dimensional representation for a ^ (3 which can 
be simply shown that it is of the form (|15p with the 
parameters given in ([7]). 

The reason that we emphasis the matrix representation 
of the algebra (fT7|) can be rewritten in the form of (fT5|) 
(which is slightly different from what was first proposed 
in (llj ) is as follows: As we have claimed in our previous 
papers, whenever the steady-state of a one-dimensional 
driven-diffusive system defined on a finite or infinite open 
lattice which evolves under the random sequential up- 
dating scheme can be written in terms of a linear su- 
perposition of Bernoulli shocks with simple random walk 
dynamics, then the algebraic relation between the oper- 
ators (when the steady-state is studied using the matrix- 
product formalism) will have a two-dimensional repre- 
sentation with a generic structure. In we have also 
proposed a general formalism by which one can simply 
find a two-dimensional representation for the quadratic 
algebra of the system in terms of the hopping rates of 
the shock front and the densities of the particles on the 
left and the right hand sides of the shock. This works if 
and only if the time evolution of the position of a prod- 
uct shock measure with a single shock front is simply a 
random walk. Moreover it has been shown, by providing 
several examples, that the conditions under which the do- 
main wall has a random walk dynamics are exactly those 
for the existence of the two-dimensional matrix represen- 
tation [8|. 

VI. CONCLUDING REMARKS 

Let us review the results of the current work. The most 
important goal in this work was answering to the ques- 
tion whether the matrix representation of the quadratic 
algebra of the system developing vis sublattice-parallel 
updating scheme has the same generic structure as we 
had proposed for the case of continuous-time updating 
scheme? By comparing our results in this paper with 
those in Q one finds that the matrix representation re- 
tains its structure even in sublattice-parallel updating 
scheme and it seems that, although we have no direct 
proof for it at the moment, the same is true for other 
updating schemes. In this direction we have considered 
a general driven-diffusive system with nearest-neighbor 
interactions which evolve vis sublattice-parallel updating 
scheme and is defined on an open lattice. If we assume 
that a single product shock measure has a simple ran- 
dom walk dynamics, generated by the transfer matrix of 
the system, then the steady-state of this system can eas- 



ily be written in terms of a linear superpositions these 
shock measures. On the other hand we have introduced 
a two-dimensional matrix representation which can gen- 
erate exactly the same steady-state. The nontrivial point 
is that this matrix representation has exactly the same 
structure that we had found for the case of continuous 
time updating scheme. 

As an evidence we have studied the TASEP under the 
sublattice-parallel updating scheme and shown that an 
uncorrelated shock can evolve in the system without re- 
quiring any constraints on the microscopic reaction rates 
i.e. the injection and the extraction rates of the par- 
ticles. The shock also reflects from the boundaries of 
the lattice with some nonzero rates. By investigating 
the time evolution equations of the shock front we have 
found that it has simple random walk dynamics. Since 
the dynamics of the shock front is quite similar to that 
of a random walker, the steady-state of the system can 
be constructed as a linear superposition of such product 
shock distributions. This could have been supposed since 
our experience with random sequential updating scheme 
had shown that in this case the quadratic algebra of the 
system should have a two-dimensional matrix represen- 
tation as it was found in As we have seen in this 
paper a two-dimensional matrix representation for the 
quadratic algebra of the TASEP under discrete time up- 
dating exists without any constraints. 
Our investigations show that the ASEP with the most 
gen eral 4-parameter open boundary conditions studied in 
[13|, can be explained using our approach provided that 
the same conditions under which the quadratic algebra of 
the system has a two-dimensional matrix representation 
are fulfilled. We have also found completely new fami- 
lies of driven-diffusive models evolving under sublattice- 
parallel updating scheme in which a product shock mea- 
sure with a single shock front has a simple random walk 
dynamics very similar to the equations ^ , provided that 
some constraints on the microscopic reaction probabili- 
ties are satisfied. We have shown that the steady-state of 
these systems can be written in terms of a combinations 
of such single shocks and at the same time the matrix 
representation of the quadratic algebras of these systems 
has the same unique structure as in ifTSj) . The details of 
these results will be presented elsewhere. 
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